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The generalised hydrodynamic theory of an electron gas, which does not rely on an assumption 
of a local equilibrium, is derived as the long-wave limit of a kinetic equation. Apart from the 
common hydrodynamics variables the theory includes the tensor fields of the higher moments of 
the distribution function. In contrast to the Bloch hydrodynamics, the theory leads to the correct 
plasmon dispersion and in the low frequency limit recovers the Navies-Stocks hydrodynamics. The 
linear approximation to the generalised hydrodynamics is closely related to the theory of highly 
viscous fluids. 

PACS numbers: 71.10.Ca, 05.20.Dd, 47.10. -hg, 62.10.-hs 
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Hydrodynamic theory of an electron gas was heuris- 
tically introduced by Bloch in 1933 |p| as an extension 
of Thomas-Fermi model. Only macroscopic variables - 
electron density n{r,t), velocity v(r, i), pressure P, and 
electrostatic potential (p{r, t), enter the theory. The set of 
equations (continuity equation, Euler and Poisson equa- 
tions) becomes complete when the equation of state is 
added, and in the original paper fit] Bloch identified P 
with the kinetic pressure of a degenerate Fermi gas. The 
Bloch's hydrodynamic theory (BHT) has been applied 
to variety of kinetic problems l|,y,§,||J6[0,@,y with mi- 
nor improvements (inclusion of exchange, correlation and 
quantum gradient corrections) ||5|,|9|]. 

From the microscopic point of view BHT cannot be 
a fully consistent theory since it extends the collision- 
dominated hydrodynamics to the electron gas where the 
collisionless (Vlasov) limit is most common. For ex- 
ample, in the plasmon dispersion law to'^ = to'^ + v^q^ 
BHT predicts for degenerate electron gas Vq = |u|. in- 
stead of a correct result |w|n [vp is the Fermi velocity) 



velocity) |17 



y,p3,h3 14 la . At arbitrary degeneracy the hy- 



drodynamics gives Wg — vj., where Vg is velocity of sound, 
whereas in the kinetic theory Wq equals to the mean 
square of the particle velocity < vl > jl^. It has been 
realized [|ri|,|2|,|3| that this discrepancy originates from 
the assumption of a local equilibrium , which underlies 
the common hydrodynamic theory [ |l6[ . The assumption 
allows to reduce the kinetic equation for the distribution 
function /p(r, t) 



dtU + -V.fp 



^-v^l^ = ^p [/p] 



(1) 



to equations for macroscopic variables n(r, t),v(r, t) and, 
in general case, temperature T{r,t). The requirement of 
the local equilibrium is fulfilled if the characteristic time 
of the process t ^ l/cj is much longer than the inverse 
collision frequency l/v and the typical length L is greater 
than the mean free path I '^ ujv {u is the average particle 



I/tv ^ ujjv < 1, ujLv <C 1. 



(2) 



In a zero oder with respect to parameters (||) Eq. (|l|) 
reduces to /p [/p] = which means that /p must have 



a locally equilibrium form, 
moving (Lagrange) frame 



The stress tensor in the co- 



^jj / .PiPjJp+mvj 

m ^ — ' 



(3) 



becomes diagonal Pij — PSij with P being the local pres- 
sure. As a result, the equations for the first three mo- 
ments of the distribution function i.e. density n, current 
j ~ nv and stress tensor Py , form the closed set of hy- 
drodynamics equations for an ideal liquid. 

Due to the high frequency of the plasma waves at least 
the first of inequalities (0) is strongly violated and the 
tensor structure of Pij as well as that of higher moments 
becomes important. In a degenerate case these tensors 
describe deformation of the Fermi sphere - the effect com- 
pletely ignored in the BHT. Within a linear response the- 
ory |1^ this effect leads to the effective shear modulus 
of a liquid. As we shall see below, this shear modulus is 
directly related to the traceless part of Pij . 

In this paper we present a generalised hydrodynamics 
which remains valid under a strong violation of both con- 
ditions (0). The theory follows from the long- wave limit 
of the Boltzmann equation and requires inclusion of the 
tensors of higher moments. We show that this leads to 
correct plasmon dispersion. 

The kinetic equation (nl) can be transformed in an in- 
finite chain of coupled equations for the moments of the 
distribution function p7[ . For the first three moments: 
density n, current j 



nv and stress tensor Pij we have 



Dtn + nVkVk = 
mnDtVi + V jPij — enVi^ = 



(4) 
(5) 



DtP^J + P,j VkVk + nk^kyj + Pkj VfeW, (6) 

+ ^kQijk — lij, 



where 



Qi]k == -^^PiPjPkfp+mv 



(7) 



is the third moment in the co-moving frame, Jy = 
■^'^pPiPjIp is the second moment of the collision in- 
tegral and Dt = dt + vV. In fl-M) the apparent conse- 
quences of the conservation of number of particles, mo- 
mentum and energy in a collision process 



^ /p = 0, ^p,;/p = 0, ;^ p2/p = 0, 



(8) 



have been used. For a charged liquid one has to add the 
Poisson equation for the scalar potential 



V (f = Anen - inpext, 



(9) 



where pext is an external charge and ip is an electrostatic 
potential. The symmetric tensor P^ (H) can be decom- 
posed into the scalar and the traceless tensor part 



P. 



P5, 



Tr-Ki 



(10) 



where P — jTrPij is the pressure of an electron gas. 

Each equation in (^^) is coupled to the next one 
through the higher moment. If conditions (g) are ful- 
filled, the higher moments are small and the infinite chain 
of equations can be decoupled (Chapman-Enskog or Grad 
methods, see p7[|). For example, as tt^ and Qijk vanish 
for locally equilibrium form of the distribution function, 
they remain small under conditions (0) . The two tensors 
describe, respectively, viscosity and thermal conductivity 
in a collision dominated liquid. However, if conditions (0) 
are violated, ttjj and Qijk as well as the higher-order mo- 
ments are not small. Yet the infinite chain of equations 
can still be decoupled if all physical quantities are slow 
varying functions of r. It is seen from Eq. (^ that the 
third moment Qijk enters only under a spatial derivative, 
hence it gives a correction proportional to 1/L. This re- 
mains true for all higher moments. Thus the truncation 
of the chain can be justified by the smallness of the gra- 
dients of the moments instead of the smallness of the 
moments itself. The dimensionless parameter of this ex- 
pansion is 



7 



Liinax{Lj, v} 



<1. 



(11) 



It is clear that both inequalities (0) can be violated while 
the condition ( pT[ ) is fulfilled. 

To show the consistency of this decoupling procedure 
we consider equations for the moments up to an arbitrary 
order in a transparent case of a one dimensional motion. 



Generalisation for a three-dimensional case is straightfor- 
ward, but calculations are lengthy and will be published 
elsewhere. 

Let the direction of the motion be along x-axis and 
hence V = xdx ■ Only diagonal x-components of all mo- 
ments contribute and the moment M^'^^ of the k-th order 
in the laboratory frame is 



mC') ^m^-''^^'"'' 



p 



Pxfp 



k>0. 



(12) 



To separate the macroscopic motion of liquid from the 
motion of particles we introduce the relative moment L^'^' 
i.e. the moment in the co-moving frame 



Lik) 



n'-'Y.p"u 



(13) 



which is related to M^*^-* as 



M«=^Ci«'L('=-'), 



(14) 



;=o 



where C^ = ^^^^^ 



^' are binomial coefficients. Apparently 

A//(o) = mn, M(i) = mnv, L^^) = p^^^ ^(3) ^ q^^^ ^nd 
L'^^ = by the definition. 

For the fc-th moment the Boltzmann equation gives 



9tM(*=) + dxM^^+^^ - k—M'^''-^^dxip = /^''\ (15) 



where I^''^ is the k-th moment of the collision integral. 
According to (|), /(o) = /(i' = 0. Substitution of (|l|) in 
(111) gives 

DfL^"^^ +L'-"^dxV^O, (16) 

L^''^DtV + dxL('^--L^"^dxV = 0, (17) 

m 

DtL'^^^ + SL^^^dxV + a,i(3) ^ j(2) ^ (18) 

DtL^''^ + {k + l)L^'''>dxV + a^iC^+i) (ig) 

Let us estimate the order of magnitude of the terms in Eq. 
(|l9|). According to Eq. ([l6|) , the spatial derivative dxV is 
of the order of the inverse characteristic time 1/t '^ uj 
hence the first two terms in Eq. dlfl) are proportional 
to uj. The right hand side of this equation is evidently 
proportional to the collision frequency i^. The last two 
terms in the left hand side of Eq. dl9) contain only spa- 
tial derivatives and are of the order of u/L. Thus contri- 
bution of L^*^' in the equation for L^^^^) contains addi- 
tional smallness 7 (pi), which means that the correction 
from L^*"'' to L^"' is of the order 7''. To obtain a theory 
valid up to 7*^ one should keep k + 1 equations for L^'^ 
(0 < Z < fc) with the contribution of the k + 1-th moment 
and the spatial derivative of L^^^ being omitted in the 



last equation. The resulting theory is a generalisation of 
hydrodynamics which is valid far from the equilibrium. 

In a second order theory one should neglect Qijk in 
Eq.(^. Using ( [lO| ) we obtain Eqs. @ and (^) in terms 
of P and TT,,; 



m nDfVi + ViP + VjVTij — enVitp — 0, (20) 

(21) 



5 2 

DtP + -PVkVk + ^TTikVkVi = 0, 



P ( V.,Uj + Vj-Wi - -% Vfc-y/c ) = /y , 



ATTij + TTijVkVk + TTjfcVfcWj + TTfcj VfcWj - -(5^ TT^fe VfcW^i 



(22) 



where Eq. (|21| ) is a trace of Eq. (|6|) and Eq. ( |22D is a differ- 
ence of ( q) a nd (pi|). The continuity equation (Q), equa- 
tions (|2C - ^2|) and Poisson equation (^ are the generalised 
hydrodynamics equations which contain two scalars n 



and P, one vector v and one traceless tensor tt, 



To 



make the system closed we also need an expression for 
lij. However, in a charged liquid contribution of lij to 
eigenmodes disappears since these modes (plasmons) lie 
in the high frequency region io/i' 3> 1 where the collision- 
less limit is reached. In this limit the linearised equations 



The first contribution in Vq equals to the square of the 
sound velocity v^ (the result of the ordinary hydrody- 
namics) and comes from the fluctuations of pressure SP. 
The second one arises from the traceless part of the stress 
tensor TTy. The sum gives a correct coefficient which in 
case of a degenerate Fermi gas reduces to the well known 
result ^Vp. It is straightforward to show that taking into 
account the third and fourth moments leads to the cor- 
rect plasmon dispersion up to q"*. 

In the low frequency (or static) limit an expression 
for collision integral is needed. We take lij in a Krook- 
Bhatnager-Gross (KBG) approximation (see, i.e. p9[|): 



/p[/p] = -^(/p(r,t)-/p^(r,i)) 



(30) 



where f^ is a local equilibrium Fermi function with 
position-dependent velocity, chemical potential /i(r, t) 
and temperature T(r,t). To satisfy the general prop- 
erties (H) the function f^ should give the same values 
of velocity v(r, t), density n(r, t) and pressure P{r, t), as 
the exact distribution function fp{r,t). 
The second moment of (pG) is equal to 



'^iP^, 



5,jPF{n,T)), 



(31) 



dtSn + noVkVk = 
mnodtVi + ViSP 
5 



VjTTy 



enoViip = 



dtSP + -PoVkVk 



dtnt 



Pn 



\/iVj 



VjV, 



(23) 
(24) 

(25) 



^S^.VkVk = 0, (26) 



where Sn and SP are deviations from the equilibrium den- 
sity no and pressure Pq. To obtain the plasmon disper- 
sion we consider a plane-wave solution e"**^'^*^''^-'. The 
velocity has only a;-component (v = xv) and only diag- 
onal component tt^x gives a contribution into Eq. (24). 
After differentiation with respect to x, Eq. ( p4| ) takes the 
form: 



^dfSn + -dlSP + -dln^x - — 5,V 
mm m 







(27) 



Combining Eq. ( p3| ) and Eqs. (^,g6|), we express 6P 
and TTxx in terms of the density fluctuations Sn: 



6P 



5/^ 
3no' 



Sn, S-Kx 



4i^ 
3no 



Sn. 



(28) 



Substitution of (|28| ) into ( |27|) together with the Poisson 
equation \^ gives the dispersion of the plasma waves: 



uj (q) =uJp+ VqQ 



where 



5Po 



4Pn 



3Pc 



3mnQ 3mno mno 



" ^<vl> 



(29) 



where Pp{n,T) is the pressure of the Fermi gas with 
the distribution function f^. By definition, this pres- 
sure equals to the exact pressure, Ppjn, T) = P = TrPij 



hence the expression in the brackets in (|3l| ) equals to nij . 
Consequently lij can be written as 



/, 



P^PF{n,T). 



(32) 
(33) 



Eqs. (|), d^jgl) with I, J from @) and P from (|3S 
constitute a closed set of equations of the generalised 
hydrodynamics in the second order of 7. Here Eq. (p3) 
does not mean the assumption of the local equilibrium, 
but simply introduces the new variable T{r, t) instead of 
P{v,t). In the low frequency limit the system is in a 
local equilibrium and r(r, t) coincides with the real local 
temperature. 

It is easy to see that Eqs. (|2(]|- p2| , p2|j3J^ ) transform into 
the common hydrodynamics in the limit lo/v <C 1. In- 
deed, in a zero order of uj/v Eq. ( p2|) leads to Tr^ = 0, 



and Eqs. (|20|) and (|2l| ) become identical respectively to 
Euler equation and to the equation for the conservation 
of energy in an ideal liquid. In the first order of to /v <^ 1 
Eq. (E2) has a solution 



= ( ViWj + VjVi - ^SijVkVk 



(34) 



which is the viscosity tensor with a viscosity coefficient 
■q — P/v. In this case Eqs.(po|, pi|) are equivalent to 
Navies-Stocks equation and the equation for the energy 
conservation in a viscous liquid. 



Though the formulas J32, |3^) were obtained using the 
KBG colhsion integral (|30|), they are more general than 
the KBG approximation itself. The reason is that in the 
low frequency range where collisions are important, the 
second moment of / is always a linear function of iiij 
| p7t . The coefficient v can thus be considered as a phe- 
nomenological parameter, which is connected to viscosity 
as 77 = P/u. 

Yet the system (0), (pQ-22) in a low frequency limit is 
not fully equivalent to classical hydrodynamics since the 
thermal conductivity contribution V'^T mm is missing 
in the energy conservation equation. Without this con- 
tribution the correct static limit Vr(r) = cannot be 
recovered. The term V^T corresponds to the correction 
of the fourth order of 7 in the continuity equation. Thus 
to include the thermal conductivity one has to consider 
the third and the fourth moments and neglect the fifth 
moment. This calculation will be presented in a more 
detailed publication. 

In the case of a dense degenerate Fermi gas T/Ep <^ 1 
[Ep being the local Fermi energy) it is unnecessary to 
consider the third and the fourth moments. At T = 



dtT^ij + tidt ( ViUj + V 



Eq. (|3^) reduces to P ■■ 
forms into condition TTijVj 



(37r2 



^/^ and (llj) trans- 
0. These two equations 
together with the system of differential equations 



'ova 



Dtn + nS/kVk = 
mnDtVi + \7iP + VjiTij — enViip = 0, 



+ P ( \/iVj + \/jVi ~ -5ij\/kVk 



provide a complete system of equations of the generalised 
hydrodynamics for a degenerate Fermi gas. 

It is interesting to note a relation of the linearised ver- 
sion of the generalised hydrodynamics in the high fre- 
quency (coUisionless) limit ( P3[ - |2q ) to the theory of elas- 
ticity. Let us introduce the displacement vector u(r, t) 
as Sn ~ — noVu. The continuity equation (|2^) gives 
the usual relation between velocity and displacement 
dfU = V. Introducing the stress tensor dij — ~6P6ij—TTij 
one can rewrite Eq. ( 24 ) as 



mnQdfUi — ^ jOi 



enQViip = 0. 



(38) 



According to the Eqs.( |25| , p6D , the relationship of the ten- 
sor ffij to Ui is the same as in the elasticity theory [Q ■ 



aij = KVkUkSij + iJi ViU 



^jUi - -SijVkUk ] , (39) 



where the bulk modulus K and the shear modulus fi of 
an electron gas are K — ^Pq and ^ = Pq. 

The inclusion of collisions violates this one-to-one cor- 
respondence to the elasticity theory. In this case the 



shear stress tensor a^ 



mined from the equation: 



■gSijakk 



SijVkUk = -VTTij, (40) 



which follows from (|2^) and (^3)- ^^ suggested in Ref. 
|18| , one may interpret this equation as if it describes the 
medium with the complex frequency-dependent elastic 
constants. However, a correspondence to so called highly 
viscous fluids is more natural. These fluids behave as 
solids at short intervals of time, but as viscous liquids on 
the large time scale [gO|. The equation for a shear stress 
tensor of such fluids has been phenomenologically intro- 
duced by Maxwell (see |20)) and exactly coincides with 
Eq. (pO[) . Hence, in a linear approximation, an electron 
gas behaves as a kind of highly viscous fluid. 
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